Effects of the crystal structure on the ferromagnetic correlations 
in ZnO with magnetic impurities 
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We study the ferromagnetism in the compound (Zn,Mn)0 within the Haldane- Anderson impurity 
model by using the quantum Monte Carlo technique and the tight-binding approximation for deter- 
mining the host band-structure and the impurity-host hybridization. This computational approach 
allows us to determine how the host crystal structure influences the impurity bound state, which 
plays an important role in the development of the ferromagnetic (FM) correlations between the 
impurities. We find that the FM correlations are strongly influenced by the crystal structure. In 
particular, in p-type (Zn,Mn)0, we observe the development of FM correlations with an extended 
range at low temperatures for wurtzite and zincblende crystal structures. However, for the rocksalt 
structure no FM correlations are observed between the impurities. In addition, in n-type ZnO with 
magnetic impurities, the impurity bound state and FM correlations are not found. 

PACS numbers: 75.50.Pp, 75.30.Hx, 75.40.Mg 
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INTRODUCTION 



Dilute magnetic semiconductors (DMS) could lead to 
new spintronic devices where both the electronic charge 
and spin can be controlled. For practical applications, 
DMS with Curie temperature, T c , above room tempera- 
ture are required. The (Ga,Mn)As is regarded as a classic 
example with robust ferromagnetism [ij , but the highest 
reported T c 's are still well below room temperature @, 
Alternatively, p-type (Zn,Mn)0 has been predicted to be 
a room temperature ferromagnet [1, Q , where ferromag- 
netism above room temperature has been observed for 
Mn doped pure ZnO 0, or Mn doped p-type ZnO 
H Ell 1 1 lL [H, HH . However, contradictory experimen- 
tal results have also been reported for (Zn,Mn)Q such as 
ferromagnetism below room temperature 141 , absence of 
ferromagnetism [l5| . spin glass behavior [l6| |. or param- 
agnetism [l7| . 

The possibility of room-temperature ferromagnetism 
motivates the theoretical research on the origin and con- 
trol of high-temperature ferromagnetism in semiconduc- 
tors, especially ZnO-based DMS, which is the subject 
of this paper. Magnetic impurities in ordinary metals 
exhibit the well-known Ruderman-Kittel-Kasuya-Yosida 
(RKKY) oscillations, which is a carrier mediated indirect 
coupling due to the Friedel oscillations of the polarized 
carriers around the impurities. When magnetic impuri- 
ties are doped into a semiconductor host, quite differ- 
ent behaviors are experimentally observed as mentioned 
above. The Haldane- Anderson impurity model had been 
introduced to study transition-metal impurities in semi- 
conductors After the discovery of DMS, the mag- 
netic properties of this model were addressed within the 
Hartree-Fock (HF) approximation , and it was shown 
that long-range ferromagnetic (FM) correlations develop 
when the Fermi level is located between the top of the 
valence band and the impurity bound state (IBS). The 
FM interaction between the impurities is mediated by 



the impurity-induced polarization of the valence electron 
spins, which exhibit an antiferromagnetic coupling to the 
impurity moments. Subsequent Quantum Monte Carlo 
(QMC) calculations [2(| on the two-impurity Haldane- 
Anderson model with the Hirsch-Fye algorithm have sup- 
ported this picture for the generation of FM correlations 
between magnetic impurities in semiconductors. Vari- 
ous other theoretical approaches have also been used to 
study magnetic impurities in semiconductors. The Zener 
model has been invoked to describe the ferromagnetism 
in semiconductors Numerical calculations based on 
local spin density approximation (LSDA) have also found 
that magnetic states and corresponding Curie tempera- 
tures in ZnO-based DMS are controlled by changing the 
carrier density or the magnetic impurity concentration 
0, SI [H, [H, HI. Within the context of DMS, the An- 
derson Hamiltonian for a semiconductor host was also 
considered by Krstajic et al. [25|, and it was shown that 
an FM interaction is generated between the impurities 
due to kinematic exchange. The role of IBS in producing 
the FM interaction in DMS was also discussed within the 
"double resonance mechanism" using HF 26]. 

In this paper, we present QMC results for the com- 
pound (Zn,Mn)0. For ZnO host, the wurtzite structure 
is the most common phase, and thus almost all experi- 
ments for ZnO are focused on this structure. But ZnO 
with the zincblende and the rocksalt structures are also 
experimentally possible in thin films and at high pres- 
sure, respectively [27j • The band structures of ZnO with 
the wurtzite, zincblende and rocksalt structures had been 
already calculated within a single set of tight-binding pa- 
rameters [27) . and we will follow these results in this 
paper to study the ferromagnetism for the compound 
(Zn,Mn)0 with experimentally determined values for the 
p-d mixing and the onsite Coulomb repulsion (28[. For 
the doped Mn 2+ impurity, we neglect the Hund coupling 
among the five occupied 3d orbitals, and, for simplic- 
ity, consider the 3c? orbitals independently. In the dilute 
impurity limit, the Haldane- Anderson impurity model is 
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invoked to describe the magnetic states of Mn 2+ ions. 
The results of the numerical calculations show that the 
crystal structure of the ZnO host strongly influences the 
energy of the IBS and the strength of the magnetic cor- 
relations between the impurities. In particular, for the 
wurtzite and zincblende structures, we find that FM cor- 
relations with an extended range develop at low temper- 
atures. However, for the rocksalt structure, FM correla- 
tions have not been observed, because, in this case, the 
IBS is located at much higher frequencies. In addition, 
only p-type (Zn,Mn)0 is found to exhibit FM correla- 
tions. 

In this paper, our purpose is to investigate the influ- 
ence of the crystal structure of the semiconductor host on 
the FM correlations between the impurities. For this pur- 
pose, we combine the tight-binding calculations for the 
host band structure and the impurity-host hybridization 
with the QMC simulations. Our impurity model is sim- 
ple, because we neglect the Hund couplings and consider 
only one of the 3c? orbitals at the impurity site. How- 
ever, this model is sufficient to demonstrate that the host 
crystal structure can be used to control ferromagnetism 
in DMS. Instead of the tight-binding approximation, the 
Local Density Approximations (LDA) can be used to cal- 
culate the host band structure and the impurity-host hy- 
bridization. In addition, it is possible to perform QMC 
simulations for all five of the impurity 3d orbitals with 
the Hund couplings. We think that this way of com- 
bining the LDA and QMC techniques can yield accurate 
predictions about ferromagnetism in DMS materials in 
the future. 



II. IMPURITY MODEL 

In order to describe the transition-metal impurities in a 
ZnO host, we use the Haldane- Anderson impurity model 
[lH which is defined by 

k,a,cr 
k,a,i,£,er 

i,£,cr i 

where cj CQ(T (ckaer) is the creation (annihilation) operator 
for a host electron with wavevector k and spin a in the 
valence (a = v) or conduction (a — c) band, and d\^ a 
(di£tr) is the creation (annihilation) operator for a local- 
ized electron at impurity site i in orbital £ and spin a with 
rii£ a = d\^ a di£ a . Here, e a (k) is the host band dispersion, 
p, the chemical potential, Vl^ka the mixing between the 
impurity and host, Ed the impurity £- level energy, and 
U the onsite Coulomb repulsion for the impurity. 

The energy bands e a (k), and the impurity-host hy- 
bridization Vi^ka will be calculated within the tight- 



binding approximation for the wurtzite, zincblende and 
rocksalt crystal structures of the ZnO host material. 
For the compound (Zn,Mn)0, the value of the onsite 
Coulomb repulsion for Mn 2+ is taken as U = 5.2eV 
by comparing with the photocmission spectroscopy mea- 
surements 28]. In addition, because the experimental 
value of Ed for Mn 2+ in ZnO host is unknown, in the 
following we use Ed = fi — U/2 so that the impurity 
sites develop large magnetic moments. The results on 
the magnetic correlations between the impurities depend 
weakly on small variations on the value of E#. 

We note that, in Eq. (1), the Hund couplings among 
the different impurity 3d orbitals is neglected. In this 
paper, we consider only one of the 3d orbitals at the im- 
purity site, since here our purpose is to demonstrate the 
effects of the host crystal structure on the FM correla- 
tions. Multi-orbital effects, where we keep all of the im- 
purity 3d orbitals and the Hund couplings, will be studied 
with QMC in a separate paper. 



III. TIGHT-BINDING APPROACH FOR THE 
ZnO BAND STRUCTURE AND THE 
IMPURITY-HOST HYBRIDIZATION 

In this section, we discuss the tight-binding calculation 
of the band-structure of the ZnO host, and the impurity- 
host hybridization. The energy bands, e a (k), of ZnO had 
been already calculated for the wurtzite, zincblende and 
rocksalt structures using a single set of sp 3 tight-binding 
parameters [27| . In this approach, the basis consists of 
one 4s and three Ap orbitals for cation Zn 2+ and three 
2p orbitals for anion O 2- . The values of the orbital en- 
ergies are E p (0) = 0.550 eV, £ s (Zn)=3.450 eV, E p (Zn) 
= 13.050 eV. In addition, the mixing values between the 
s, p orbitals of Zn 2+ and the p orbitals of O 2- are taken 
to be (spa) = 2.965 eV, (ppa) = 4.324 eV, and (ppn) = 
-1.157 eV. 

Using these tight-binding parameters and keeping all 
of the branches within the sp 3 basis, we have reproduced 
the band structure of ZnO. In Figs. QJa)-(c), we have 
plotted the branches near the semiconductor gap. Here, 
we observe that, for the wurtzite and zincblende struc- 
tures, the top of the valence band is located at the T point 
with a direct gap of 3.45 eV. For the rocksalt case, the 
top of the valence band, located at the L point, is at 0.6 
eV, while the T point is at -0.53 eV. Hence, for the rock- 
salt structure the system has an indirect semiconductor 
gap. 

Next, we discuss the calculation of the impurity-host 
hybridization within the tight-binding approximation. 
Once a substitutional impurity Mn 2+ is introduced and 
takes the position of a Zn 2+ cation, the 3d orbital £ of 
Mn 2+ will mix with the neighboring 2p orbitals of O 2- . 
The mixing matrix element V^ka = ( ( P^(^)\H\'^ a (k)) has 
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FIG. 1: Energy bands of the ZnO host near the semiconductor 
gap with (a) wurtzite, (b) zincblende, and (c) rocksalt crystal 
structures obtained using the tight-binding approximation. 
These results were reproduced using the parameters given in 
Ref. IH. 



the form of 



(2) 



where cp^(i) is the impurity 3d-state at site i, and 5 , Q (k) is 
the host state with wavevector k and band index a, which 
is expanded by atomic orbitals y p (n) with orbital index 
p and site index n. Here, N is the total number of host 
lattice sites, and a ap (\i) is an expansion coefficient. For 
the p-d mixing integrals of (ip^(i)\H\ip p (ri)) , p represents 
the three 2p orbitals of 2 ~, and £ denotes the five 3d or- 
bitals of Mn 2+ . As shown by Slater and Koster [2!| , these 
fifteen mixing integrals can be expressed by only two in- 
tegrals (pda) and (pdn) and direction cosines I, m and n 
in the two-center approximation. In this approach, the 
integrals (pda) and (pdn) are taken as fitting parameters 
which can then be determined by other exact results or 
by comparison with the experimental data. For the com- 
pound (Zn,Mn)0, the value of (pda) = -1.6 eV has been 



estimated by comparin g w ith the photoemission spec- 
troscopy measurements [28[, while the value of (pdw) is 
always determined by the relation (pdn) = — (pda)/2.16 
[30l ] . We will use these values in the rest of this paper. 

It is established that, for the wurtzite or zincblende 
crystal structures, the p-d mixing is dominated by the 
occupied t2 g (xy, yz, zx) orbitals because of the tetrahe- 
dral crystal field, while for the rocksalt structure, the p-d 
mixing comes mainly from the occupied e g (x 2 — y 2 , z 2 ) 
orbitals due to the octahedral crystal field [31(. In addi- 
tion, the Hund coupling among the 3d orbitals of Mn 2+ is 
neglected in this paper. In addition, for simplicity, here 
we consider only one of the t2 g orbitals (£=xy orbital 
here) for the (Zn,Mn)0 with wurtzite and zincblende 
structures, and only one of the e g orbitals (£ = x 2 — y 2 
orbital here) for the (Zn,Mn)0 with rocksalt structure. 
Hence, we treat the 3d orbitals of Mn 2+ independently. 
In this paper, we are mainly interested in the effects of 
the host crystal structure, and we will treat the multi- 
orbital effects in a separate paper. 

Figure[5]displays results on the the p-d mixing function 
V^(k) defined by 



^(k)=(X)|^a(k)| ! 



iV2 



(3) 



where only one of the Mn 2+ 3c? orbitals, labeled by £, 
is considered. In Eq. ([3|), the summation over a is per- 
formed only over the valence bands (Fig. 2(a) and Fig. 
2(c)) or the conduction bands (Fig. 2(b) and Fig. 2(d)). 
Here, V^(k) is plotted along various cuts in the Brillouin 
zone for wurtzite, zincblende and rocksalt crystal struc- 
tures. Figures 2(a) and (b) show V for a Mn 3d(xy) 
orbital when the ZnO has the wurtzite structure. Here, 
we observe that, at the T point, the total hybridization 
of the xy orbital with the valence bands is about three 
times larger than that with the conduction bands. For 
the wurtzite and the zincblende structures, the semicon- 
ductor gap edges are located at the T point, hence the 
value of V near T will be particularly important in de- 
termining the energy of the IBS and the strength of the 
magnetic correlations between the impurities. Figures 
2(c) and (d) show that, for the case of a Mn 3d(xy) or- 
bital in ZnO with the zincblende structure, the total hy- 
bridization with the valence bands is also stronger than 
that with the conduction bands near the T point. Fig- 
ures 2(c) and (d) also show results for a Mn 3d(x 2 — y 2 ) 
orbital in the rocksalt structure. Here, we see that the 
total hybridization with the valence and the conduction 
bands vanish at the T point. However, at the L point, 
where the top of the valence band is located, the total 
hybridization with the valence band is the stronger than 
in the rest of the cases. 

The results on the frequency of the IBS and the 
strength of the FM correlations depend sensitively on 
the value of hybridization near the gap edges, which we 
will discuss in the next section. We note that, in turn, 
the hybridization depends strongly on the values of the 
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FIG. 2: Hybridization function V^(k) of a single Mn 2+ or- 
bital with the valence bands ((a) and (c)) or the conduction 
bands ((b) and (d)). Here, results are shown for a Mn 2+ 
3d(xy) orbital in the wurtzite or zincblende ZnO and for a 
Mn 2+ 3d(x 2 — y 2 ) orbital in the rocksalt ZnO. In obtaining 
these results, we have used mixing parameters with values 
taken from Ref. [281 . 



mixing parameters {pda) and (pdn) of the tight-binding 
approach. 



IV. QUANTUM MONTE CARLO RESULTS 

In this section, we present results on the impurity mag- 
netic correlations, which were obtained using the Hirsch- 
Fye QMC technique [33 ]. The input parameters for the 
QMC simulations were calculated with the tight-binding 
approach described above. The following results were 
obtained with more than 10 5 Monte Carlo sweeps and 
Matsubara time step At = 0.225. 

We first discuss the local moment formation at an im- 
purity orbital. For this purpose, we have performed QMC 
simulations to calculate ((M z ) 2 ), where 



M z 



(4) 



is the magnetization operator for a single 3d(£) orbital at 
the impurity site i. We have performed the calculations 
for a single 3d orbital because in this paper we are mainly 
interested in the effects of the host crystal and we neglect 
the multi-orbital effects, which we will treat in a separate 
paper. Hence, ((M z ) 2 ) represents the square of the local- 
moment for a single 3c? orbital added to the ZnO host. 

Figure 3(a) and (b) shows results on ((M z ) 2 ) for a 
id{xy) orbital in wurtzite and zincblende structures and 
a 3d(x 2 — y 2 ) orbital in the rocksalt ZnO. In Fig. 3(a), 



((M z ) 2 ) versus the chemical potential p is plotted for 
< /i < 0.35eV and in Fig. 3(b) for < p < 4eV. These 
results are for temperature T = 100A". In these figures, 
we observe discontinuities at 0.12eV, 0.20eV and 1.6eV 
for the zincblende, wurtzite and rocksalt structures, re- 
spectively. 

According to the Hartee-Fock and QMC calculations, 
the presence of a discontinuity in ((M z ) 2 ) versus p, im- 
plies the existence of an IBS at this energy [l9|, [2(| • A 
step discontinuity develops in the magnitude of the mag- 
netic moment as p increases through loibSi because im- 
purity spectral weight is induced in the semiconductor 
gap at to ibs for sufficiently strong hybridization. When 
the IBS is occupied, the spin polarization of the host 
split-off state, which is due to the impurity-host mix- 
ing and at the same energy as the IBS, will cancel the 
spin polarization of the valence band. Thus antiferro- 
magnetic couplings between the polarized host carriers 
and the impurities disappear. This causes the FM in- 
teraction between the impurities, which is mediated by 
the polarized carriers around the magnetic impurities, to 
vanish. In this paper, we are studying how this mech- 
anism of ferromagnetic correlations is influenced by the 
crystal structure of the the ZnO host. The variation in 
the values of ujibs for different crystal structures seen 
in Fig. 3 is clearly a consequence of the differences in 
the energy bands and the impurity-host mixing. In addi- 
tion, the values of the hybridization with the bottom of 
the conduction bands are weaker, and hence we do not 
observe bound states near the bottom of the conduction 
band. 

The value of loibs plays an important role in determin- 
ing the strength of the FM correlations which develop be- 
tween the impurities when the IBS is unoccupied. Within 
the Hartree-Fock approximation and for a semiconductor 
host with constant density of states po and semi-infinite 
bands [n| , the range of the ferromagnetic (FM) correla- 
tions between the impurities is given by 



1 



(5) 



when the IBS is unoccupied (0 < p < to ibs)- However, 
when the IBS becomes occupied (luibs < m); the FM 
correlations become weaker. The QMC calculations per- 
formed for a two-dimensional (2D) semiconductor host 
with quadratic quasiparticle dispersion confirm this pic- 
ture [2fJ|. These QMC calculations show that, in 2D, 
coibs increases as the strength of the hybridization grows. 
The maximum range of the FM correlations decreases as 
coibs increases in agreement with the HF calculations. 
However, in the three-dimensional case, the IBS does not 
exist, if the hybridization is smaller than a critical value 

In the following calculations, we will see that the mag- 
netic correlations between the impurities are sensitive to 
the location of the chemical potential with respect to the 
IBS. An IBS with shallower position (cof BS ~ 0.1 eV) 
is obtained for (Zn,Mn)0 with the zinblende structure, 
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while a much deeper IBS (uifgg ~ 1.6 eV) is found for 
the rocksalt case. Based on the previous HF and QMC 
calculations [III H3] , it is reasonable to expect that the 
FM correlation range for (Zn,Mn)0 with the wurtzite 
structure will be shorter than that with the zincblende 
structure, and will be much longer than that of the rock- 
salt case. 

In addition, it should be pointed out that the IBS 
for (Zn,Mn)0 with wurtzite and zincblende structures 
only exists near the top of valence band, while no IBS is 
found near the bottom of conduction band. If, instead 
of Mn 2+ , another transition metal ion TM 2+ is substi- 
tuted into the ZnO host, then the only different quanti- 
ties would be the p-d mixing parameters {pda) and (pdn) 
within this framework. In particular, for various hosts 
and transition-metal impurities, the p-d mixing param- 
eters, which are consistent with the experimental mea- 
surements, are mostly in the range —1.6 eV< (pda) <- 
1.0 eV (1, H, while (pdn) = -(pda)/ 2. 16 We 
have checked the IBS for (Zn,TM)0 for various p-d mix- 
ing values in the above mentioned range. We find that, 
when smaller values are used for the p-d mixing, the IBS 
shifts towards the top of valence band, while at the same 
time no IBS develops near the bottom of the conduction 
band. Since the doping of the transition metal ion TM 2+ 
into the ZnO host does not itself introduce carriers, the 
hole and electron carriers in (Zn,TM)0 are associated 
with additional the acceptor or donor defects, respec- 
tively. This implies that for p-type (Zn,TM)0, the FM 
correlations will develop when < \i < ujibs, while for 
the n-type (Zn,TM)0, no FM correlations are expected 
due to the absence of the IBS. In fact, the existence of 
p-type rather than n-type (Zn,Mn)0 with wurtzite struc- 
ture has recently been pointed out based on the analysis 
of the experimental measurements [l2l [l3j . 

We next display results on the impurity-impurity mag- 
netic correlation function (Mf M|) versus the impurity 
separation R/a in Fig. 0J These results are for two Mn 2+ 
3d orbitals, which are of the same type, added to the ZnO 
host. Clearly, in a more realistic calculation of FM cor- 
relations between the impurities, it would be necessary 
to consider the magnetic correlations between different 
types of 3d orbitals. However, our purpose in this paper 
is to investigate the role of the host electronic structure. 
In Fig. 4, the temperature T = 200K, R = \R r - R 2 \, 
and a is the lattice constant. In our chosen directions, 
the distance between the lst-nearest-neighboring Mn 2+ 
is a for the wurtzite structure, while it is a y/2/2 for the 
zincblende and rocksalt structures. As shown in Fig. Ufa) 
for the wurtzite structure, the impurity spins exhibit FM 
correlations at /i = 0.0, but it is short range. By in- 
creasing /j, to 0.15 eV, the range of the FM correlations 
becomes longer extending to the 3rd-nearest neighbor- 
ing Mn 2+ . Further increasing /i, the FM correlations be- 
comes weaker. This is because the IBS of (Zn,Mn)0 with 
wurtzite structure becomes occupied as \i changes above 
0.15 eV, as seen in Fig. [3] For the zincblende structure, 
similar results are obtained as shown in Fig. HJb): the 
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FIG. 3: Square of the magnetic moment at the impurity site, 
((M z ) 2 }, versus the chemical potential fi for different energy 
intervals at T = WOK. These results were obtained for a 
single Mn 2+ 3d(£) orbital added to the ZnO host. For the 
wurtzite and zincblende crystal structures, we have considered 
a £ = xy orbital, while for the rocksalt case an x 2 — y 2 orbital. 
The vertical dashed line denotes the bottom of the conduction 
band. 



impurity spins exhibit FM correlations extending to the 
lst-nearest neighboring Mn 2+ at /x = 0.0. Increasing fi 
to 0.1 eV, the FM correlations extend to the 4th-nearest 
neighboring Mn 2+ . Further increasing /j, to above 0.1 eV, 
the FM correlations become weaker. This is because the 
IBS of (Zn,Mn)0 with the zincblende structure becomes 
occupied as \i > 0.1 eV, as displayed in Fig. [3] We 
note that the compound (Zn,Mn)0 with the zincblende 
structure seems to possess a longer range for the FM 
correlations than that with the wurtzite structure, while 
almost all existing experiments for ZnO are focused on 
the wurtzite structure. Not surprisingly, we have not 
observed magnetic correlations for the rocksalt structure 
in Fig. HJc). This is because the position of IBS for 
(Zn,Mn)0 with rocksalt structure is too deep as seen in 
Fig. [31 thus the FM correlation range 0, [2(| is shorter 
than the first nearest-neighbor distance. 

We next study the temperature dependence of the fer- 
romagnetic correlations between the impurities in differ- 
ent crystal structures. Figure 5 displays the T depen- 
dence of (Mf M| ) for impurities which are second-nearest 
neighbors in wurtzite, zincblende and rocksalt crystal 
structures. For wurtzite and zincblende structures, these 
calculations have been performed for the 3d(xy) orbitals, 
and for the rocksalt case they are for the 3<i(x 2 — y 2 ) or- 
bital. In addition, for the wurtzite and zincblende struc- 
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FIG. 4: Impurity-impurity magnetic correlation function 
(MfMf) versus distance R/a for (a) the wurtzite, (b) 
zincblende, and (c) rocksalt structures. Here, each impurity 
site consists of a single Mn 2+ 3d orbital. For the wurtzite 
and zincblende crystal structure, we have considered a 3d(xy) 
orbital at the impurity sites, while for the rocksalt case a 
3d(x 2 — y 2 ) orbital. 



tures, the chemical potential is taken as /i = 0.15 eV 
and 0.1 eV, respectively, where the IBS is unoccupied 
and the longest range of the FM correlations occurs as 
shown in Fig. [U The results for the rocksalt case are for 
/i = 1.5eV, so that the IBS is unoccupied here also. We 
observe that, as T decreases from 800-ftT down to 200K, 
the FM correlations grow for the wurtzite and zincblende 
structures. In this T range, the FM do not develop for 
the rocksalt case. 



V. DISCUSSION 

For the ZnO host with the rocksalt structure, it is use- 
ful to mention previous studies [U [H, HI] which are 
related to the indirect energy gaps obtained with the 
sp 3 tight-binding calculation shown in Fig. [ljc). Us- 
ing the Hartree-Fock method, Jaffe et al. [35| pointed 



A 



V 



0.2 



0.1 



0.0 



Between the 2nd -nearest Mn 




200 



400 

T[K] 



600 



800 



FIG. 5: Impurity-impurity magnetic correlation function 
(MfMf) for second nearest-neighbor impurity sites versus 
temperature T. Here, for the wurtzite and zincblende crystal 
structures, we have considered 3d(xy) orbitals, while for the 
rocksalt case 3d(x 2 — y 2 ) orbitals. 



out that for the rocksalt ZnO, the point symmetry in 
rocksalt structure does not allow for mixing between the 
Zn 3d orbitals and the O 2p orbitals at the T point, but 
mixing is allowed elsewhere in the Brillouin zone. Thus 
the valence band maximum shifts away form the T point, 
so that the gap becomes indirect. As mentioned by Skin- 
ner and LaFemina in Ref.[27j], where we follow their sp 3 
tight-binding calculation for the ZnO host in present pa- 
per, this effect also exists even when the sp 3 tight-binding 
model does not explicitly include the Zn 3d orbitals, be- 
cause the Zn 3d character has been implicitly included 
in the hopping integrals through the tight-binding inter- 
polation of the bulk ab initio pseudopotential bands [36| 
that Skinner and LaFemina used to derive them. 

However, there also exist calculations for the rocksalt 
ZnO, which are based on the quasiparticle approach (37j 
and the pseudopotential method [38| . yielding a direct 
energy gap at the T point. In order to understand the 
results observed for (Zn,Mn)0 with the rocksalt struc- 
ture, we have considered the case of MgO, whose crys- 
tal structure is also rocksalt but has a direct energy gap 
at the r point 39]. Using the tight-binding parameters 
given in Ref. [39| , where the host basis consists of one s 
orbital for Mg 2+ and three 2p orbitals for O 2- , we have 
calculated the p-d mixing for the compound (Mg,Mn)0 
as in the case for (Zn,Mn)0. We found that, near the 
r point, the p-d mixing for (Mg,Mn)0 with the rocksalt 
structure is also close to zero, where d stands for the five 
3d orbitals. These results imply that the host compound 
with the rocksalt structure does not mix with the doped 
3d orbitals at the T point regardless of whether there is 
a direct or indirect energy gap. In fact, this is gener- 
ally true if the host with the rocksalt structure can be 
described using the sp 3 tight-binding approach. At the 
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r point, because of the orthogonality of the s and p or- 
bitals, we find that, for a fixed band a, the coefficients 
a ap in Eq.((2]) have only one non-zero value, for example, 
a. Thus the p-d mixing at the T point has the form 

V a (0) = a£> f (i)|/%p(n)) = - n) (6) 

n n 

where £ represents one of the five 3d orbitals, p denotes 
one of the s, p x , p y , and p z orbitals. In addition, here, 
the summation is taken over the nearest neighbors, and 
E( n) is given using the notation of Slater and Koster 
[29( | . which can be expressed in terms of the non-zero 
(pda) and (pdir), and the direction cosines I, m, and n. 
For the rocksalt structure, it can be shown that the value 
of Vq)(0), Eq. ©, vanishes for any £ and p orbitals, when 
the summation is performed over the six nearest neigh- 
bors. 

It is obvious that the magnetic correlations between 
the impurities depend on the position of the IBS, while 
the IBS is closely related the host band structure as well 
as the host-impurity mixing parameters. As mentioned 
above, for the rocksalt ZnO doped with magnetic impu- 
rities, the IBS will shift to the top of valence band if we 
decrease the values of the p-d mixing parameters (pda) 
and (pdir). Thus the detailed band structure calcula- 
tions, such as LDA, are needed in the future to combine 
with QMC so that the ZnO-based DMS can be described 
more accurately. 

VI. SUMMARY AND CONCLUSIONS 

In summary, we have studied the ferromagnetism for 
the compound (Zn,Mn)0 with different crystal struc- 
tures in the dilute impurity limit based on the Haldane- 
Anderson impurity model. The band structures of the 
ZnO host were calculated using the sp 3 tight-binding pa- 
rameters from Ref. [27l |. and the p-d mixing parame- 
ters and the onsite Coulomb repulsion U were obtained 
from comparisons with the photoemission measurements 
on (Zn,Mn)0 [Hj]. The QMC calculations show that the 
magnetic correlations between Mn 2+ impurities in ZnO 
are strongly affected by the host crystal structure. For 
the wurtzite and zincblende structures, the ferromagnetic 



correlations are found, and their range extends up to the 
third or the fourth nearest neighbor sites at low tem- 
peratures. On the other hand, for the rocksalt struc- 
ture, no magnetism has been found even between the 
nearest-neighbor impurities. In addition, only p-type 
ZnO doped with magnetic impurities is found to have 
ferromagnetism. 

In this paper, our main purpose has been to investi- 
gate the effects of the crystal structure on the IBS and 
the FM correlations for the ZnO host. For simplicity, we 
have neglected the Hund coupling among the Mn 2+ 3d 
orbitals. In addition, we have considered only one of the 
five 3d orbitals at the impurity sites. This is clearly an 
over simplification, however, it does allow us to demon- 
strate the role of the host crystal structure in determin- 
ing the energy of the IBS and the maximum range of 
the FM correlations. Currently, we are in the process 
of performing QMC calculations which include all five 
of the 3d orbitals at the impurity site without neglect- 
ing the Hund couplings. In addition, we are using the 
LDA technique, instead of the tight-binding approxima- 
tion, to calculate more precisely the host band structure 
and the impurity-host hybridization. We think that this 
LDA+QMC approach can be used to make accurate pre- 
dictions about the FM correlations for different sets of 
DMS materials. 

It is generally agreed up on that the ferromagnetism in 
the DMS can be controlled by changing the type of the 
transition-metal impurities or the host semiconductor as 
well as the occupation of the IBS. The results presented 
in this paper suggest that the host crystal structure can 
also be used in the search for high-T c ferromagnetism in 
DMS. 
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